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Abstract 

We consider zero-temperature transitions from conformal to non-conformal phases in quantum 
theories. We argue that there are three generic mechanisms for the loss of conformality in any 
number of dimensions: (i) fixed point goes to zero coupling, (ii) fixed point runs off to infinite 
coupling, or (iii) an IR fixed point annihilates with a UV fixed point and they both disappear 
into the complex plane. We give both relativistic and non-relativistic examples of the last case in 
various dimensions and show that the critical behavior of the mass gap behaves similarly to the 
correlation length in the finite temperature Berezinskii-Kosterlitz-Thouless (BKT) phase transition 
in two dimensions, £ ~ exp(c/|T — Tcj 1 / 2 ). We speculate that the chiral phase transition in QCD 
at large number of fermion flavors belongs to this universality class, and attempt to identify the 
UV fixed point that annihilates with the Banks-Zaks fixed point at the lower end of the conformal 
window. 
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I. INTRODUCTION 



The renormalization group (RG) underlies our understanding of second-order phase tran- 
sitions, with critical points being identified with fixed points of the appropriate RG equa- 
tion pp. Near the phase transition the characteristic energy or momentum scale m (the 
inverse correlation length) goes to zero as m ~ \ot — a*\ v , where a is a parameter that can 
vary continuously, and a = a* is the location of the critical point. 

In this paper, we argue that there is wide class of phase transitions in which the correlation 
length behaves very differently, vanishing exponentially on one side of the phase transition, 
while being strictly zero on the other side 



m ~ A uv 6(a* — a) exp , c > . (1) 

This peculiar behavior — where all derivatives of the correlation length with respect to a 
vanish at the critical point — has been observed before in the Berezinskii-Kosterlitz-Thouless 
(BKT) phase transition in two dimensions [2]; therefore we will refer to eq. ([TJ as "BKT 
scaling." The BKT transition is a classical phase transition in two dimensions that can be 
described in terms of vortex condensation. It arises due to the competition between the 
entropy of a single vortex and the binding energy of a pair of vortices, both which scale as 
log-R, R being the size of the system. While this transition is peculiar to two dimensions, 
we will show that the mechanism underlying BKT scaling from an RG point of view is far 
more general, and is one of three generic behaviors that can occur when a system in any 
dimension makes a transition from a conformal to a non-conformal phase. In particular, as 
we will show, it follows when an IR fixed point of the system merges with a UV fixed point. 
In this language it is easy to see why BKT scaling can be found in a wide variety of systems. 
The basic mechanism can be illustrated with a simple model with a dimensionless coupling 



g depending on an external parameter a, for which the /5-function takes the form (Fig. 1(a) ) 



dg 

P(9;a) = = (a -a*) - (g-g*) 2 , (2) 

where t — In//, fi being the renormalization scale. For (a — a*) > 0, the fixed points for this 
system (zeros of 0) are given by 

g± = g* ± y/a-a* , (3) 

where g_, g + correspond to IR and UV fixed points respectively, each describing a conformal 
phase of the theory 1 . As a decreases, these two fixed points approach each other until they 
merge at g± = g* for a = a*. For a < a* the solutions to (3 — are complex, and the theory 
no longer has a conformal phase. 



1 By IR and UV fixed points we mean zeros of the /^-function which are attractive or repulsive in the IR 
respectively. 
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(a) (b) 

FIG. 1: (a) A toy /3-function. For a > a* there are fixed points at g± which are UV- and IR-stable 
respectively; these fixed points merge at g* for a = a*, and disappear for a < a*; (b) The RG flow 
of the coupling g as a function of t = In fi in the non-conformal phase, with (iuv - ^IR.) °t 1 /V a * ~ a - 



To see that fixed point merger generically gives rise to BKT scaling, consider the case 
where a is slightly below a*, and that at a UV scale Auv the coupling takes an initial value 
9uv < 9*- On scaling to the IR, the coupling then flows to larger values, lingering near 
9 = 9* where the /3-f unction is small, and then blowing up quickly, denning an intrinsic 
IR scale Air, which is insensitive to the initial value guv- This behavior is displayed in 



Fig. 1(b) The scale Air will characterize the longest correlation lengths in this theory, and 



can be computed by integrating eq. 



exp [t m - tuvj = exp 



A 



uv 



gm dg 



where we have assumed |^ir,uv ~ 9*\ ^ \ a ~ a *\ ■ 

How general is this mechanism of fixed point annihilation? Suppose a system has a 
nontrivial IR fixed point at g — g(a) whose location depends continuously on a parameter 
a, and that at a critical value a = a* there is a phase transition where conformality is lost. 
At this phase transition, the /^-function must somehow lose a zero. This can come about in 
three ways: 

(A) g can decrease until it merges with the trivial fixed point at g = 0, giving rise to a 
trivial, asymptotically unfree theory; 

(B) g can run off to infinite coupling and disappear; 

(C) g can merge with a UV fixed point, as in our toy model, giving rise to BKT scaling. 

Examples of scenarios (A) and (B) are afforded by supersymmetric QCD (SQCD). At large 
number of colors N c , the parameter x = Nf/N c may be treated as continuous, where Nf is 
the number of quark flavors. It has been shown by Seiberg [31 Sj that SQCD is conformal 
in the window 3/2 < x < 3. For x just below 3, the theory has a Banks-Zaks fixed point 
at weak coupling [5]; approaching x = 3 from below, this fixed point merges with the trivial 
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fixed point at g = 0, and for x > 3 the theory is in the asymptotically unfree "free electric 
phase." This is an example of mechanism "A" above. In contrast, at the lower end of the 
conformal window at x — 3/2, SQCD goes from a strongly coupled conformal theory when 
x > 3/2 to a "free magnetic phase" when x < 3/2. In the free magnetic phase, the Coulomb 
force between charges takes the form e 2 ln(Ar)/r 2 where A is associated with the Landau 
pole of the dual magnetic theory. The log behavior of the coupling can be explained by a 
/5-function which is negative and approaches zero as (3 ~ — 1/g for large g. Thus it appears 
that conformality in the electric description is lost via mechanism (B). [Yet, since in the dual 
magnetic theory conformality is lost via mechanism (A), it would appear that scenarios (A) 
and (B) can describe the same physics in terms of different degrees of freedom.] 

In this paper we give several examples of theories which exhibit the mechanism (C) of fixed 
point merger and BKT scaling. Following our RG analysis of the original BKT transition, 
we analyze the quantum mechanical example of a 1/r 2 potential in d dimensions, which can 
be solved nonperturbatively and which exhibits the phenomenon of fixed point merger. We 
show how this analysis has many parallels in the AdS/CFT correspondence [6j [7J E], and 
that loss of conformality via fixed point merger is analogous (if not holographically dual) to 
the instability of AdS space at the Breitenlohner-Freedman (BF) bound [9]. 

Our next example is a relativistic theory of gauged fermions confined to a defect. Here 
a perturbative analysis near d = 2 dimensions reveals fixed point merger and BKT scaling. 
A rainbow approximation to the gap equation gives qualitatively similar results. 

One of the motivations for this paper is to understand the chiral phase transition that 
happens in (nonsupersymmetric) large-A^ QCD when the number of flavors of massless 
fermions Nf varies. As with SQCD, we know there exists a conformal window for QCD 
in the parameter x = Nf/N c where the upper end occurs at x* = 11/2, near which the 
Banks-Zaks calculation is perturbative and reliable. For decreasing x conformality must 
eventually be lost, since for small x chiral symmetry breaking is expected. We speculate 
that the phase transition at this lower boundary of the conformal window occurs due to 
fixed point merger. This suggestion is not new: it has been advocated before by Gies and 
Jaeckel based on the results from the functional RG approach [10]. If this picture is correct, 
then near the transition the chiral condensate must exhibit BKT scaling. Incidentally, this 
exponential behavior is also typically found when one solves the gap equation obtained from 
(an unsystematic) truncation of the Schwinger-Dyson hierarchy (see, e.g., [TT t IT2 T IT3 T [T4"l ITS] . 
and [16J for further references). 2 A priori, the relationship between the RG picture of merging 
fixed points and the gap equation is not obvious; however our analysis of relativistic defect 
fermions yields the same result in the regime where both approaches are reliable. 

If QCD does indeed exhibit BKT scaling, then our arguments suggest that within the 
conformal window there exists another theory, QCD*, which is defined at the UV fixed point. 
We conclude with speculations about this theory. 



2 In this context, BKT scaling is sometimes called "Miransky scaling. 
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II. THE BKT PHASE TRANSITION 



The BKT phase transition [T?l IT8] is due to the deconfinement of vortices in the XY 
model at a critical temperature T c , above which the theory is conformal. The behavior 
of the correlation length eq. below the phase transition can be understood from the 
appropriate RG equation [IS]. We can exploit the equivalence between the XY model and 
the zero temperature sine-Gordon model in 1 + 1 dimensions: 

L = |(^0) 2 -2^cos0, (5) 

where T corresponds to the temperature of the XY model in units of the spin coupling. Near 
the phase transition, it is useful to use the variables u = 1 — 1/ 8nT and v = 2z/TA 2 — where 
A is the UV cutoff associated with the vortex core — in terms of which the perturbative 
/3-functions are 

(3 U = -2v 2 , p v = -2uv . (6) 

Changing variables to v + u = r and v — u = 2w, one sees that rw invariant under RG flow, 
and the running of r is governed by 

dr 

/3(r; wt) = n— = —2wr — r 2 . (7) 
d\x 

This /^-function has exactly the quadratic form of our toy model eq. with the substitution 

(a - a*) -> -2wt, (g - g+) -> r. (8) 

However, this /3-function is only valid for small r and w, so the region about r = is 
excluded for fixed wt, as shown in Fig. [2j Because of the excluded region, the physics for 
the BKT model is slightly different than for the toy model: in the non-conformal phase 
(wt > 0), instead of starting from the left of r* = in the UV and flowing to the right in 
the IR, the system starts at the top of the hill just to the right of r* and flows to the right 
in the IR. While it may appear that this requires a fine-tuned initial condition for r, that 
is not the case in terms of the u and v variables. Starting the flow near r = gives a factor 
of 1/2 in the exponent for the correlation length relative to the expression eq. Q: 



;bkt 



A ~ e 7T/{2V ^ ) . (9) 



The critical temperature is found by solving wt = 0; expanding about T = T c yields the 
familiar BKT result 

£bktA ~ e «l\T-^< 2 , (10) 
where b' is a nonuniversal number that can be expressed in terms of z and A. 
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FIG. 2: The function /3(r) in the vicinity of r = for the BKT transition eq. ([7|; the gray region 
is outside the realm of validity of the calculation. 

III. A NONRELATIVISTIC EXAMPLE: QUANTUM MECHANICS IN 1/r 2 PO- 
TENTIAL 



It is well known that the solutions for a quantum particle in a potential 

V(r) = a/r 2 (11) 

possess conformal symmetry when the potential is repulsive or weakly attractive (a > a*), 
but that for sufficiently attractive potential (a < a*), conformality is lost and the poten- 
tial has discrete bound states. 3 For a range of a, the zero-energy, s-wave solution to the 
Schrodinger equation for two particles with mass m — 1 in d dimensions interacting via the 
potential V(r) is given by 

V = c_r"- + c + r*+ , z/± = -i!_^ ±v /^^; ) a ^JA^L. (12) 
This solution is valid for a in the range 

a* < a < a* + 1 ; (13) 

for a < a* the above solution becomes complex, and the Hamiltonian does not have ground 
state, while for a > a* + 1 then z/_ < — d/2 and the r u ~ solution is not normalizable near 



r ~ 0. Within the range eq. (13), if either c + or c_ vanish, then the solution is scale- 
invariant. Solutions for which both c + and c_ are nonzero define an intrinsic length scale, 
L = (c + / 'c_) 1 /( I/ - _1/ +) and therefore do not exhibit conformal invariance; however in this case 
the solution always approaches c + r u+ for large r (since v + > vJ) and so we can identify 
the c_ = solution with an IR attractive fixed point and the c + = solution with a UV 
fixed point, in a manner we can make precise. Arranging to have one of these solutions or 
the other requires different boundary conditions at the origin, so we see that the theory is 



3 There is a vast literature on the 1/r 2 potential. For textbook treatment, see Ref. [2D]; for an early 
reference, see [2T]; for relatively recent RG treatments see [531 122 [25] and references therein. 
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actually not well denned with the potential eq. (11), but that it must be augmented by a 5 
function at the origin which controls the boundary condition at r = 0: 



V(r) = a/r 2 - g5 d (r) . (14) 

We will show that the coupling g obeys an RG equation analogous to our toy model eq. (J2J), 
and that the two conformal solutions c_ = and c + = will correspond to two different 
fixed points of the coupling g. As a approaches a* from above, we will show that the two 
fixed points merge, g± — g*, at a value for g* which we will compute. For a < a* a UV cutoff 
must be imposed on the theory in order to have a ground state and an IR scale emerges 
which is related to the UV cutoff through the BKT scaling formula eq. Q. We show this in 
two different ways: first we perform a nonperturbative analysis, and then we use Feynman 
diagrams in a perturbative calculation in 2 + e dimensions. Both calculations shed light on 
the relativistic example we provide later, and on our conjecture about the behavior of QCD 
as a function of the number of flavors. 



A. A nonperturbative calculation 



1. The exact wavefunction and energy 

To solve the Schrodinger equation exactly for two particle scattering via a 1/r 2 potential 
in d dimensions we need to regulate the singularity at r = 0. We choose to do so by 
considering the potential 



V(r) 



a/r 2 



r > r , 



(15) 



where r n 1 will serve as the cutoff A 



uv- 



g/r 2 , r < r , 

At low energy, there is a region tq < r <C 1/ \fE where the Eip term in the Schrodinger 

(16) 



equation can be neglected, and for a > a* we find the solution eq. (12) 

ifj = c_r u ~ + c + r u+ , 
with the ratio c + / c_ given in terms of Bessel functions as 



c_ 



-r, 



1 = 



Ji=l (y/g) 



(17) 



The quantity (c + /c_) is a dimensionful quantity characterizing this solution; by requiring 
that it does not change as we change the UV cutoff r , we arrive at the exact /3-function for 
7 (defining RG time t = — lnr ): 



$7 
di 



-(7 + f+)(7 + = (« - «*) - (7 - l*f 



(18) 



7 



d-2\ 2 d-2 d-2 



with 

a* = - [~^~) ' 7 * = ' 7± = ~ v t = ± ■ (19) 

We recognize this to be the same /^-function as our toy model eq. ^ with fixed points at 



7±; referring to eq. (17) we see that the IR fixed point corresponds to 7 = 7_ and c_ = 0, 
while the UV fixed point is associated with 7 = 7+ and c+ = 0. 

For general d and a < a*, scaling solutions do not exist; physical quantities, such as 
the bound state energy, depend on the UV cutoff. Motivated by the discussion of coupling 
constant flow in our toy model, we know that physical quantities will be insensitive to the 



value of 7 (the UV coupling) so long as 7 < 7*, as seen in Fig. 1(b) So we take 7 — > —00, 
which is reached in the limit of a hardcore repulsive potential for r < r , g — > —00. The 
ground state wavefunction is then described by the Bessel function ip(r) = r~( d ~ 2 ^ 2 K ir) {kr) 
with i] = \Ja* — a and the boundary condition ip(ro) = 0. For small real rj we can solve for 
k and find the binding energy 

9 1 ( 2tt 



£ = ^= exp +0(1) , (20) 

Note that this scale for the binding energy is easily attained from the RG analysis as 

B ~ Aj R = fl e n^/^ 2 = JL e -2-M^ . (21) 

If one takes 7 to be arbitrarily close to 7* in the UV, and then takes a — > a* then the 
binding energy goes to zero, but the exponent is only half as large, 

B~\exp(- * ). (22) 

recalling the result for the BKT transition. 



2. Onset of the Efimov effect 



Although the /3-function in eq. (18) takes the same form as the toy /3-function in eq. ([2]), 
eq. (17) implies that the coupling g is a multi-valued function of 7, with I7I — > 00 identified 



with the zeros of J d-2 (y/g) ■ Therefore for a < a* our RG equation actually describes limit 
cycle behavior: as 7 runs from —00 to +00 in RG period T = — f dj/fiy = n/y/a* — a, 
g runs from one Bessel function zero to the next. It follows that there is not just one IR 



scale defined by this RG flow, as in Fig. 1(b), but an infinite number of such scales, each 
successively smaller than the previous by a factor of exp[— n/^a* — a}. This behavior can 
explain the Efimov effect in 3-body bound states. 

The classic Efimov effect [27] concerns the system of three identical bosons. When the 



scattering length between two bosons becomes large, the three-body system develops a 
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series of ever shallower bound states. This occurs because the three particles interact via an 
a/r 2 potential for r 3> Tq, where r$ is the two-body effective range, and a a fixed number 
satisfying a < a*. These systems require a 3-body interaction, and the renormalization of 
this interaction exhibits the limit-cycle behavior discussed above. The infinite tower of IR 
scales is associated with the infinite number of "Efimov states" below threshold, exhibiting 
a geometric spectrum [26J. Such states have been observed in systems of trapped atoms 
tuned to a Feshbach resonance. 

Three degenerate bosons tuned to infinite scattering length (so-called "unitary bosons") 
do not have a variable a parameter; in order to see a transition very similar to what happens 
at a = a* we need a case when the Efimov effect appears as one changes a tunable parameter. 
This is realized by nonrelativistic fermions at unitarity with different masses for two spin 
components, M (heavy) and m (light). The Efimov effect occurs in the p-wave channel for 
two heavy and one light fermions if M/m > 13.6 [28J. 

It is known that for 8.6 < M/m < 13.6 one can additionally fine tune the three-body 
interaction to resonance [29]. From our point of view, the two theories with and without 
fine-tuning in the three-body channel correspond to the UV and IR fixed points. When 
M/m — ► 13.6, the two fixed points approach each other: the difference between theories with 
and without 3-body fine-tuning becomes smaller and smaller. Finally when M/m > 13.6 the 
fixed point completely disappear, and an energy scale appears in the problem: the ground 
state energy of the three-body bound state. 



3. Operator anomalous dimensions at the IR and UV fixed points. 

We can gain insight about the two fixed points by looking at the dimension of the op- 
erators. Let us consider the two-particle operator ipip. According to the operator/state 
correspondence developed in Ref. [30J, one can find dimensions of this operator by putting 
two particles in a harmonic potential. The Hamiltonian of the system is given by 

H = - \v\ + V(\ ri - r 2 |) + \uj\t\ + r\) . (23) 

In terms of the center of mass coordinate R and relative coordinate r, the Hamiltonian can 
be rewritten as H = Hr + H r where the ground state energy of Hr equals dui/2 and 

H r = -V 2 r + V{r) + ^V, (24) 



where the potential is given in eq. (15). The ground state wavefunctions and energies for 



this Hamiltonian for g tuned to one of the fixed points g± is easily seen to equal 

i>± = e-^l'r^ , Ef =(± + „ ± ) uj (25) 

in the limit r — > 0. (Recall that the fixed points g± correspond to solutions ip = r u± in the 
absence of the harmonic potential). Therefore the total ground state energy is 

E ± = (d + v±) uj , (26) 
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and so the scaling dimensions of the two-particle operator tptp are 

A± = (d + v±) = V^^, (27) 

where A + and A_ are the operator dimension at the IR and UV fixed points respectively. 
We emphasize the fact that 

(A+ + A_) = d + 2 (28) 



in any spatial dimension d and for any in the range a* < a < (a* + 1) in eq. (13), with 
A + = A_ = (c?+2)/2 at a = a*. Note that (d+2) is the scaling dimension of a nonrelativistic 
Lagrange density, since time has twice the scaling dimension as space; we return to this 
below, when we discuss the AdS/CFT correspondence. 

B. The renormalization group: e expansion 

When we consider relativistic quantum field theories a nonperturbative solution will not 
be available, and we will have to rely on either perturbation theory, or a truncation of the 
Schwinger-Dyson equations. It is therefore instructive to examine a perturbative analysis 
of the 1/r 2 potential. We have seen that the doubly degenerate fixed point at the phase 
transition occurs at coupling 7 = 7* = (d — 2)/2; we therefore start with the action in 
d = 2 + e spatial dimensions, where perturbation theory can correctly describe the phase 
transition. In order to facilitate the use of Feynman diagrams, we write the theory in second 
quantized form with a contact interaction, 



S = jdt rf d x (i^dfip - 



- j dtd d xd d y^(t,x)^(t,y) 



; i(>(t,y)1>(t,x), (29) 
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i x ~y| 

where the factor of tt in the contact interaction is chosen for future convenience. The 
Feynman rules are as follows: 



Propagator 



to - p 2 /2 



(30) 



Contact vertex 



"Meson exchange" 



i-KQH \ (31) 



— A- (32) 
e |q| e 
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FIG. 3: Two diagrams contributing to the (3- function in Eq. (33). Note that the second diagram 
is a tree diagram. 



Note the unusual 1/e pole in the "meson propagator." It arises because the Fourier transform 
of a 1/r 2 interaction is log divergent in d = 2 dimensions. 

It is easy to see that a does not get renormalized (as one would expect, being the strength 
of a nonlocal interaction); however the coupling g runs. From the above Feynman rules, the 
/^-function for g arises from the sum of the tree graph and the one-loop graph shown in 
Fig. |3j with the result 



I3(g; a) = ^ = eg-^ + 2a = 2^a + jj--(g-eY, (33) 

which we recognize to be equivalent to our toy model, up to an unimportant rescaling of g 
by 2, with 

g* = e, a* = - — , (34) 
Note that our perturbative expansion is justified for small e, but a* coincides with the exact 



result in Eq. (19). For a > a* the /3-function has two zeros: g± = g* ± 2y/a — a*. At 
a = 0, g_ = is the IR stable fixed point, corresponding to a noninteracting theory — for 
the generic short-ranged potential, low-energy scattering is trivial; g + = 2e corresponds to a 
fine-tuned potential with a bound state at threshold (e.g., an infinite scattering length). As 
one decreases a the two fixed points approach each other, merging at g± = g* when a = a*. 

For a < a* the potential requires a cutoff and has a bound state; we can estimate the 
size of the bound state to be given by the correlation length £ = in eq. Q; this gives a 
binding energy B ~ Aj R , or 

B~A 2 VY exp(- J^— ) ■ (35) 

Note that this formula is independent of e and therefore appears to be independent of 
dimension. In fact the above estimate is verified in the nonperturbative calculation of the 
previous section. 

An unusual feature of our calculation of the /3-function (Fig. [3]) is the contribution from a 
tree graph. We close this section by noting that a more conventional calculation is obtained 
by making the following change of variable: 

9 = ~9-^. (36) 
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FIG. 4: Diagrams contributing to the perturbative /3-function for g. 

Then the RG equation becomes 

dq ~ q 2 n ~fOL\ 1 / 2a\ 2 . . 

i = ^-^ + 2 9 (-)--( KT y (37) 

The first term on the the right hand side comes from the engineering dimension. The other 
terms come from the diagrams as in Fig. |4| All diagrams now have loops. This is a more 
natural approach from the point of view of Wilsonian RG, where one looks at the logarithm 
in the momentum integral instead of the 1/e poles. But we emphasize that the two RG 
equations lead to the same physical consequences. 

1. Summary of the QM example 

Before proceeding, we summarize our findings. 

• The theory has two fixed points, IR and UV, when a > a* = —(d — 2) 2 /4. 

• When there are two fixed points, the dimensions of the scalar operators at the IR and 
UV fixed points are A + and A_, and they satisfy A + + A_ = d + 2. 

• When a < —1/4, the fixed points do not exist and the theory develops a bound state 
energy which scales as Ay V exp(— 2n/^a^ — a). 

IV. A HOLOGRAPHIC PERSPECTIVE 

As far as we know, there is no simple holographic dual description of quantum mechanics 
with 1/r 2 potentials, nor of the field theoretical models considered later in this paper. How- 
ever, holography provides an interpretation of conformality loss which turns out to be very 
useful in developing our intuition about such phase transitions: the loss of conformality can 
be associated with the violation of the Breitenlohner-Freedman (BF) bound. 

A. The conformal phase: pair of theories 

In our RG discussion, for a > a* there are two CFTs that merge into one at a = a*. 
This situation is reminiscent of what occurs in holography [51] : a higher dimensional theory 
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containing a scalar field with mass m 2 in the interval — d 2 /4 < m 2 < — d 2 /4 + l corresponds 
to two different boundary theories in which the dimensions of the operator O dual to have 
two different values 

A ± = ^± 1 -Vd 2 + Am 2 = ^±u. (38) 

The main point here is that in the asymptotics of the scalar field near the boundary z = 
of the AdS space, <f>(z) = c_z A ~ + c + z A+ , one can interpret c_ as the source coupled to O, 
and c_|_ as its expectation value, and vice versa. 

Instead of repeating the discussion in Ref. [31], we illustrate its main points in a simple 
model. In this model, one sees that the theory with [O] = A_ can be obtained from the 
theory with [O] = A + by adding to the Lagrangian a term O 2 with a fine-tuned coefficient 
(in other words, we will go "against the RG flow," cf. Ref. |32j where one follows the RG 
flow from the UV fixed point to the IR fixed point). 

Consider a massive scalar field in AdS^+i space. We use Euclidean signature in this 
subsection, so the metric is 

ds 2 = -^r(dz 2 + dx^dx^ 1 ) . (39) 
z 2 

We will set the radius of the AdS space R — 1. The action for the scalar field = <ft(z,x) is 

S = i jdzd d x y/g^d^ut + m 2 <f) 2 ) - [d d x J(x)0(e, x) 

2 J e + J (40) 

= l - Jdzd d x^ l [z 2 (d z 0) 2 + z 2 (d^) 2 + m 2 <j ) 2 ]--^ Jd d xJ{x)<f>{e,x). 

In our model, this action is taken as the the definition of the CFT. This CFT "lives" on the 
boundary in the sense the external source J couples only to the field at some small z — e, 
with 1/e playing the role of the momentum UV cutoff. The operator 0(x) that J couples to 
is defined as 0(x) = e~ A+ 0(e, x). The extra power of e is chosen so that subsequent results 
have a regular e — > limit. 

We assume a large N parameter so that one can use the saddle point approximation, in 
which satisfies the field equation 

0" _ ^1^' _ g 2 _ + eA-ij^z - e ) = o , (41) 

z z z 

where we have changed to momentum space. We assume qe 1, i.e., q is much smaller 
than the UV cutoff. To completely specify the solution we impose two boundary conditions. 
Near z = there are two possible solutions to this equation, ~ z A± where A± are defined 
in Eq. (38). We require that 

<f) = c z A +, z^O, (42) 

i.e., we require to follow the most regular asymptotic behavior at small z. We leave the 
boundary condition at z — > oo for later discussion. Equation (42) is valid for z < e, but due 
to the insertion of a source at z = e, contains both asymptotics once z is larger than e, 

= c+ z A+ +c_z A ", eKz^q' 1 . (43) 
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Clearly, c_ is proportional to the source J. Matching boundary conditions one finds 

J 



C- 



(44) 



From the point of view of the interior region z > e, Eqs. (43) and (44) effectively fix 
the boundary condition near z — e. Here we obtain a key ingredient of the AdS/CFT 
prescription: the coefficient in front of the z A ~ part of the field is the source coupled the 
the operator 0. 

The coefficients Cq in Eq. (42) and c + in Eq. (43) can be determined only after the 
boundary condition at z = oo is fixed. We can relate the expectation value of O with c + : 

4>(e) J 



(A- 



A _VA+-A_ 



(45) 



Therefore, up to a singular contribution, the expectation value of O is related to c+. The 
two-point function (00) is then 



(00) 







(0)1 



J dc+ 
+ 



(46) 



dJ""" J 2ve 2u ' dJ 

Let us now impose the boundary condition at z — » oo. To ensure finiteness of the action, 
it is sufficient to require 

<f>(z,x) — >0, z — ► oo . (47) 
The saddle point solution is now completely determined 



' De- u K v (qe)z d / 2 I v (qz), z < e 
De- v I u (qe)z d ' 2 K v (qz), z > e 



(48) 



with D = Je " '. This solution corresponds to 

J 

c_ = — , c + = 
2u + 



m-u) 



(49) 



z/ 2 r(z/)2 1 + 2i - 

The two-point function (00) is proportional to q 2u , consistent with dimension of O being 
A + =d/2 + v. 

Now we turn on a deformation O 2 with a coefficient that will be fine-tuned to get another 
conformal field theory. The action is now 



S 



dz d d x —— [z 2 {d z <j)) 2 + z 2 (d^f + m 2 2 ] - / d d x 



X 

2^' 



(e) + 



Let us first set J = 0. The field equation is 



+ 



d 



+ 



m 



A 



5{z-e)<f) = 0. 



(50) 



(51) 



z z* e 

One can integrate this equation from z = to larger z. For z < e, <p is purely z A+ , and 
for z > e it becomes a mixture of z A+ and z A ~ , the relative weight of which depends on A. 
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The most interesting value of A is when is purely z A ~ for z > e. This happens when A is 
fine-tuned to the critical value 

A = A + -A_. (52) 

There is a quantum-mechanical interpretation of this fine-tuning. If one identifies z as 
the radial coordinate r of a two-dimensional space, then Eq. (51) is the radial Schrodinger 
equation for the wave function ip = z~ d ^ 2 (f) of a particle moving in a potential which is a 
sum of a 1/r 2 piece and a delta-shell piece, 

V(r) = ^ - -S(r - e) . (53) 
r z e 



with — q 2 playing the role of the energy. The value (52) corresponds to the case when the 
potential has a zero-energy bound state. 

Let A be fine-tuned to this value, and turn on the source J. Using the asymptotics 
<fi ~ z A+ for z < e and integrating the field equations passed z = e, we find that for z > e, 
the coefficient c + is now proportional to J: 

The expectation value for O is now related to c_, 

(0)j = c_ . (55) 

The assignment of source and expectation value is reverse to the case A = 0. If one imposes 
the boundary condition <p(z) —* when z — > oo, then the solution to Eq. (51) is given by 
Eq. ((481), but now 

Je u 

D = l-2uI v (qe)K u (qe) ' (56) 

The solution corresponds to 

J 2 2 - 1 i» J 

C+ ~~2~v> C ~ " r(l - u) ■ { } 

In particular (OO) ~ q~ 2u , corresponding to [O] = A_ = d/2 — v. 

Thus, in this simple holographic model, the UV stable fixed point of the CFT with 
the fine-tuned O 2 interaction corresponds to the same bulk theory, but with the opposite 
assignment for the source and the expectation value. 



B. Below the Breitenlohner-Freedman bound 

Here we speculate on the fate of the bulk theory with a scalar with m 2 below the 
Breitenlohner-Freedman (BF) bound —d 2 /4. The most interesting case is when m 2 is only 
slightly below the BF bound, where the boundary theory is approximately conformal over 
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a large energy range. The dual bulk description should involve a spacetime that is approxi- 
mately AdS, cutoff both at the UV and the IR by the respective "walls." 

First, for the set up with a scalar m 2 below the BF bound, there must be an UV cutoff in 
the theory. For example, the theory with a dual description can arise as a low-energy limit 
of another theory whose UV is free of any instability. Let us model that by imposing a hard 
cutoff on the AdS space, and impose a boundary condition on the scalar at the cutoff. 
The precise form of the boundary condition is not important, for definiteness we take it to 
be Dirichlet: 4>(z\jy) = 0. 

One expect that a IR scale will be generated by the condensation of 0. We model that 
scale very roughly by another, IR, cutoff at Zir, and impose another Dirichlet boundary 
condition there. 

Now let us look at the field equation for 0, 

0" - —0' - ^0 - g 2 = . (58) 

z z 2 

Changing variables to = z^ 1 ^ 2 ^, this equation becomes 

- V + — —i' = • (59) 

This equation, with the boundary condition on at Zir and zjjy, gives us an infinite tower of 
particles. The mass square of the particles in this tower is the eigenstate of a particle in a 
one-dimensional potential which is a/r 2 enclosed between two infinite walls at z^y and zir. 
The condition of absence of tachyon is equivalent to the condition that the potential does 
not contains a negative-energy eigenstate. This requires the interval between the two cutoff 
is not too large, 

ln^< - I =. (60) 

^UV V m BF ~ m2 

In a more realistic setup where the scale zj^ appears dynamically, one can expect that it 
appears at the scale requires for preventing a tachyon, 

zm ~ zuv exp ( J = ] . (61) 



\Anfp - m- 



We expect that this situation is rather generic in holography. It would be interesting to 
construct an explicit solution in string theory which exhibits the BKT scaling. 

V. A RELATIVISTIC EXAMPLE: DEFECT QFT 

In this section we consider a relativistic quantum field theory that exhibits the phe- 
nomenon of fixed point annihilation. The example resembles QCD with large number of 
flavors, but the phase transition occurs in the regime of weak coupling. 
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We consider a theory of a fermion living on a d-dimensional membrane, and interacting 
through a SU(iV c ) gauge field that lives in (3+1) dimensions. We shall assume that the 
SU(iV c ) gauge coupling does not run; it can be easily accomplished by taking the gauge field 
to be part of a conformal field theory, say, the M = 4 super- Yang-Mills theory. The most 
interesting case is d = 3 (=2 + 1), which was analyzed by Rey using the Schwinger-Dyson 
approximation and gauge/gravity duality [33J. Here we shall take d = 2 + e (i.e., (1 + e) + 1) 
to take advantage of a small parameter e < 1. 

The action is 



S 



d d x {ijyfW + ^VM*) - 4 / d * x 4 



(62) 



and we assume there is a UV cutoff A. The <i-dimensional photon propagator is obtained 
by integrating out the 4d propagator over transverse directions, 

d 2 ~ e q± —i 



(2tt) 2 - q 2 



9,u-(l-0 



q 2 



where q 2 = q 2 — q\. For small e, the result is 

%9p.v I 1 



1 



(63) 



(64) 



2vre \(-q 2 y/ 2 

Note that the dependence on the gauge parameter £ disappears in the small e limit (if one 
assumes £ ~ 1). 



A. Schwinger-Dyson treatment (rainbow approximation) 



Before going to the RG treatment, we review how the chiral phase transition is found 
using the gap equation. This treatment is very similar to that used in QCD [Til E]- The 
lowest-order gap equation is (Fig. [5J 

d d q 



iE(p) 



(27T) 



q 2 — Xr(g) + le 



(65) 



Inserting the photon propagator (63 ) and performing a Wick rotation, the equation becomes, 
for small e, 



E(p) 



9 2 C A 
47r 3 e 



d d q 



E(g) 



|p — g| e g 2 + S 2 (g) 



^ 2 



i 



2AT, 



(66) 



P- 3 




FIG. 5: The one-loop graph that contributes to the gap equation 
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where the integral is taken in Euclidean space. It will become clear later that the dominant 
contribution to the integral comes from the regions p <C q and p ^> q, with p ~ q giving a 
subleading contribution. Changing variables to 



x = In 



P 



y = hi 



Q 



m m 

where m = S(0) will be the mass gap, the Schwinger-Dyson equation becomes 



Ux) 



g 2 c A 

27r 2 e 



dy [ e - €(x - y) - e-< Xm -^\L{y) + dy[l- e-< Xm -^]E(y) 



UQ 



(67) 



(6* 



where x m = ln(A/ra) and A is the UV cutoff. Differentiating Eq. (68) over x, we find 

9 2 C A 



S'(x) 



2tt 2 



dye-« x -^(y) 



E"(x) = f X dye-< x -y^(y) - ^E< 



x) 



(69) 
(70) 



from which we find that S satisfies the differential equation 

Z"(x) + eX'(x) + ^Z(x) = 

with boundary conditions 

S'(0)=0, Z(x m ) = 0. 
The solution to the equation is 

E(x) =m e-^ CQS( ";^ 
coso 

with 



K 



9 2 C A 



(71) 



(72) 



(73) 



(74) 



2tt 2 4 ' 

When e and k are small, £ varies slowly on the logarithmic scale, which validates the 
assumption that the integral in Eq. (66) is dominated by regions where p and q are very 
different. 

The boundary conditions imply 



tan 8 



2k 



cos(nx m — 8) — 



from which one finds 



x , ■ 



n H — + arctan 



2k 



(75) 



(76) 
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where n is an integer. The solution with n = corresponds to the biggest gap and is favored 
energetically. The dynamically generated mass gap is 



m ~ A exp 



1 /7T e 
— — h arctan — 
k \2 2k 



K 



9 2 C A 
2tt 2 



(77) 



So we find that there is a phase transition occurring at 



gl 



2C A ' 



(78) 



and the critical behavior of the gap near g = g* conforms with BKT scaling. 



B. RG treatment: beyond the rainbow 



The RG equation can be written in a way very similar to the RG equation for the 
QM example with 1/r 2 potential. One introduces an extra four-fermi interaction into the 
Lagrangian 



S 



C , T 



d a x i^dud + ipi^A^ - -{^nyy - 



d A x F^F" + ■ 



(79) 



The tree level one-gluon exchange contains a 1/e factor from the gluon propagator (64) and 
contributes to the beta function for c: 



f3( c ) = ec - —c 2 



Nc 
2tt 



3_ 
2tt 



The phase transition occurs at g = g* where /3(c) has a double zero, 



vr 2 e 2 



N c 



When g > g*, we need to solve the RG equation, 

dc 



(3(c) 



(80) 



(81) 



(82) 



with the boundary condition that the bare four-fermi coupling is zero at the UV cutoff, 
g(A) = 0. The solution is 



ne 2ti 



-re tan 



k In 5 



K 



g 2 N c 

An 2 



N N 

where 5 = arctan(e/2«;). The coupling constant becomes infinite at 

1 /7T 



(83) 



m = A exp 



1 /7T 

-A-2 +6 



(84) 
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We find that in the limit N c — > oo, the result from the RG approach coincide with what 
is obtained from the Schwinger-Dyson approach. However, for finite N c the results of the 
two approaches are different. This is not unexpected, since the RG sums up a wider class 
of diagrams than the gap equation. 

For g < (?*, there are two zeros of the /^-functions 

c± = y (vre =f y/n^ - g*N c ) . (85) 

On the other hand, the scaling dimension of operator ip-ip is 

N e i 

A±[^] = 1 + e - -^c± = 1 + - ± yJ^-g>N c . (86) 

We find that 

A + + A_=2 + e = d. (87) 
up to possible corrections of order e 2 . 



C. Summary of the relativistic example 



The lessons we learn from the relativistic examples are very similar to the nonrelativistic 
example: 

• g < g*: two fixed points. The dimensions of the operator ipifj at the two fixed points 
satisfy A + + A_ = d. 

• g > g*' no fixed points, gap formation, BKT scaling. 



VI. QCD AT LARGE N c AND N f . 



We now turn our attention to the most interesting, and most difficult, example; the 
chiral phase transition in QCD with large iV c and Nf. Denote x = Nf/N c . We consider the 
Veneziano limit iV c — > oo, Nf — > oo, x fixed, and denote the (rescaled) 't Hooft coupling as 

a ~ 92Nc 

as ~ Wf ■ (88) 

The beta function of QCD in this regime is 

P(a*) = ~\ [(11 - 2x ) a l + ( 34 - 13a a ' + •••]• (89) 

For x < 11/2 the theory is asymptotically free. If x is slightly below 11/2 the beta function 
has a nontrivial zero which is still in the perturbative regime: 

a s * = ^-(ll-2x). (90) 
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This is the Banks-Zaks (BZ) fixed point [5]. In the IR, the theory is an interacting CFT. 
This fixed point moves to strong coupling as one makes 11/2 — x ~ 1. For small i,a;<l, we 
believe that the theory has chiral symmetry breaking and a confinement scale. It is natural 
to assume that there is a critical Nf/N c ratio a; cr i t at which the chiral condensate goes to 



zero. 



If the picture emerging from the previous examples also holds for QCD, then conformality 
is lost when the BZ fixed point annihilates with another UV fixed point. Therefore, we 
predict that when x is slightly larger than x cr i t , QCD has an UV fixed point, in addition to 
the IR fixed point (and the free UV fixed point). 



Am 



\\)\\) 2- 




Free fermions 



QCD 



Xbz 
= 11/2 



FIG. 6: A possible picture for the QCD chiral phase transition in Nf/N c . The lines denotes the 
dependence of the dimension of the chiral condensate ipip at the fixed point. The solid line is the 
IR fixed point, and the dashed line is the UV fixed point. 

The situation is illustrated in Fig. [6j The UV fixed point called QCD*, is a different CFT 
compared to the usual IR fixed point; for example, the dimension of the operator iptp should 
be different between the two fixed points. 

What is the nature of QCD*? It could be that the /3-function for the QCD gauge coupling 
simply has an unstable fixed point at strong coupling. However, this picture implicitly 
assumes that the set of relevant operators in QCD* consist of just kinetic terms for the gauge 
fields and fermions, as is the case at weak coupling. At strong coupling other operators could 
be relevant as well, and guided by our defect QFT example of Sec. [Vj it is natural to consider 
the possibility that a chirally symmetric four-fermion operator is relevant in QCD* 



C 



QCD* 



C 



QCD 



(91) 



and that the unstable fixed point exists at some value {g*, C*} in the two-dimensional space 
of couplings. By analogy with Sec. |Vl we then expect that the beta function for c contains 



4 One assumes that some UV scale, e.g., the scale Aqcd defined from one-loop running, is fixed when x is 
changed. 
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linear, quadratic and constant terms, 

(3(c) = 7l c - l2 N c c 2 + 7o3 2 , (92) 

where the linear 71c term is due to the anomalous dimension of the four-fermi operator, the 
quadratic c 2 term is due to the one-loop graph involving two four-fermi vertices, and the 
constant g 2 term is due to, e.g., one-gluon exchange graph. 5 This is essentially the picture 
advocated by Gies and Jaeckel [ID]- 6 The constants 70, 71, 72 depend on x = Nf/N c ; and 
x cr it is where /3(c) has a double zero. 

We do not know where in x the fixed point QCD* exists. A particularly interesting possi- 
bility is that QCD* exists at weak coupling, say near x = 11/2. As described subsequently, 
we find many theories similar to QCD* in the perturbative regime, but none of them possess 
the full chiral symmetry of QCD, and hence cannot be QCD*. The possible "phase diagram" 
is illustrated in Fig. [6| the line corresponding to QCD* does not continue to the vicinity of 
x = 11/2. 

In the rest of this section, we shall be looking for perturbative UV fixed points that flow 
to the BZ fixed point. Models A and C below have been considered in Ref. [36J. 



A. Model A 

As the first step, we try to find a fixed point when one operator, ipif), has dimension 
different from its dimension at the BZ fixed point. From AdS/CFT experience, we expect 
that operator dimensions at the two fixed points satisfy 

A+ + A_ = d = 4. (93) 

At the BZ fixed point A + « 3, therefore at the QCD* fixed point, A_ rs 1, which means that 
this operator is almost a free scalar. This suggests that we look for QCD* in the following 
Lagrangian, which will be called model A, 

£modcl A = £ Q CD + \{d^f ~ ~ ^ . (94) 



It is convenient to define the rescaled couplings, 

_ g 2 N c _ y 2 N c N f ~ _ XN c N f 



(4tt) 2 ' ay (4tt) 2 ' " (4tt) 2 



(95) 



5 The constant term can be of a different power of g, but it does not affect the argument. 

6 Beta functions containing three terms of the almost the same origin arise in orbifolds of Af — 4 super- 
Yang-Mills theory for the coefficients of double-traced operators [34] . The constant term also appears in 
the running for the Landau liquid parameters in the RG treatment of color superconductivity |35j . 
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These constants will be found to be 0(N®) at the fixed point. In this regime, the beta 
functions are 

!3 aa = -\ [(11 - 2x)a 2 s + (34 - 13x)a 3 s ] , (96) 

Pa y = -6a s a y + 2a 2 , (97) 

fa = -12a 2 , + 4%A (98) 

The fixed point for the gauge coupling, a s *, is the same as at the BZ fixed point to leading 
order in 1/N C . The fixed point for the Yukawa and four-scalar couplings are 

StJtj;); j A 3(2^ Qft^ . (^99) 

Thus, model A has a perturbative fixed point. 

We can compute the dimension of the scalar operators at the fixed point, 



jbz = 3 - 3a s *, (100a) 
A [0] mo d e i a = 1 + a y * • (100b) 

Notice that at in the model A fixed point, the operator ipip is replaced by the operator 0. 
We see here that 

A[#] BZ + A[0] modcl A = 4 , (101) 



which coincides with our expectation (93). From the QM intuition, we may expect that 
when A + = A_ = 2, the BZ fixed point and model A become identical. 

According to our expectation, the new fixed point should be an UV fixed point, and 
that there exist a deformation of this fixed point that leads to the BZ fixed point. The 
deformation is provided by the scalar mass term m 2 2 . This deformation is relevant if A[0] 
is less than 2. If such perturbation is present, the scalar decouples below a certain energy 
scale, leaving the theory to be in the BZ fixed point, as expected. 

It might seem that the way model-A Lagrangian was introduced, with an extra scalar 



field and Yukawa interaction, is very different from the way done in Eq. (91 ). It seems that if 
we want to change the dimension of ipip, then one should introduce a four-fermi interaction 
into the QCD Lagrangian: 

Anodcl A = ^QCD + C(^V) 2 . (102) 

We argue here, however, that the two forms of the Lagrangian are just two different rep- 



resentations of the same fixed point; with Eq. (94) being the weak-coupling representation 
near the upper end of the conformal window (x = 11/2), and Eq. (102) being the more 
useful representation near the lower end (x = x C rit)- Indeed, the propagator of at the 
IR fixed point is g 2A ^ _4 where A<^ is the dimension of 0; near the lower end of the confor- 
mal window A = 2, therefore the scalar propagator is almost momentum independent, and 
the scalar-mediated interaction between fermions becomes point-like. This is similar to the 
equivalence between Nambu-Jona-Lasinio and Yukawa models in dimensions between 2 and 
4 [33138]. 
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B. Model B 

Model B is similar to model A, except there are now two scalar fields, 

£> = £qcd + i(^) 2 + \(d,<p 2 ) 2 - JLjfa + ^ 7 5 2 )^ - ^(4>1 + <Pl? ■ (103) 

The Lagrangian preserves vector SXJ(Nf) and axial U(1)a (more precisely, the nonanomalous 
discrete subgroup of it). The behavior of model B is exactly like in model A: there is a fixed 
point for y and A; and the running of g is not altered in large N c , Nf regime. 



C. Model C 

Both model A and B preserves only a small subset of the SU(A^)xSU(iV/) chiral sym- 
metry of QCD. The simplest Lagrangian which preserves chiral symmetry is 

C = C Q cb - y(ipt A tlj(f) A + i$t A j 5 ipir A ) + Trd»&d fl $ - Ai(Tr<&t<&) 2 - A 2 Tr($+$) 2 , (104) 

where $ = (0 A + iic A )t A , A = 0, . . . Nf are flavor Gell-Mann matrices, normalized so that 
Tr(t A t B ) = \5 AB . We need to find the fixed point of this theory. This fixed point should 
have only one IR unstable direction corresponding to the mass term for 

The RG equations for g and y can be read out from Refs. (39J, SOj E]. We need the 
two-loop beta function for the gauge coupling (as the one-loop contribution has a small 
coefficient near x = 11/2), but for the Yukawa and scalar couplings one- loop beta function 
suffices. Moreover, the two-loop beta function for g and the one-loop beta function for y 
does not contain scalar self-couplings, thus one can first solve for the fixed point for g and 
y, and then look for the fixed points of scalar couplings. In this model we define a y as 



y 2 N c 



«v = 773! ■ ( 105 ) 



The beta functions are 



2,. „ x o 2 



(3 as = - 2x)a 2 s - -(34 - 13x)a 3 s - 2x 2 a 2 s a y , (106) 

(3 ay = -6a s a y + 2(1 + x)a 2 y . (107) 

We work around x — 11/2. At fixed a s , there is a zero of f3 y , but there is no fixed point of 
both beta functions. Therefore model C does not have a perturbative fixed point. 



D. Model D 

The difference between model C and models A, B is that we introduce 0(N 2 ) scalars into 
model C, while there are only 0(1) scalars in models A, B. As a result, the beta function 
for the gauge coupling changes, and there is no longer a fixed point. 
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Our last model, model D, interpolates between models B and C. We introduce couplings 
to 2M 2 scalars that preserve a SU(M)xSU(M)xSU(fc) subgroup of the chiral symmetry 
group, with M = Nf/k, 



£ = £qcd - yip° l t A /3 ((j) A + i7 5 7r A )^f + scalar terms, 



(108) 



where a, (3 runs 1 . . . M, i runs 1 . . . k, A runs 1 . . . M 2 . Model B corresponds to k — Nf 



and model C to k — 1. We redefine a„ to be 



y 2 kN c 
(4tt) 2 



The beta functions are now 

/?«. = -2a 2 s 
(3 ay = 2a 



ll-2a: 34-13x 
^ + ^ 



X 

0,8 + T2 a y 



-3a s + 11 + 



x 

V 2 



For x slightly below 11/2, the model has a fixed point for any integer k larger than 1, 



2k 2 



11 11 -2x 



2k 2 



[11 - 2x) 



(109) 

(110) 
(111) 

(112) 



25A; 2 -44 3 ' v 25k 2 - 44 

but there is no perturbative fixed point or k = 1 (model C). (The fixed point values for the 
four-scalar couplings can also be found.) 

Therefore, there exist theories that preserve part of the chiral symmetry of QCD and 
flow to QCD by a relevant deformation, but we have not succeeded in finding a theory that 
plays the role of QCD* in the perturbative regime. This does not mean QCD* does not 
exist; in fact we will give arguments, largely based on holography, that it does exist in the 
nonperturbative region. 



Operator dimensions in model D. 

The dimension of the scalar operators (p A , n A in model D is 

A[0] | model D = 1 + tiy* | model D • (H3) 



On the other hand the dimension of tptp at the BZ fixed point is given in Eq. (100a). Taking 
the sum of the dimensions, we find 

A[^]| BZ + A[0]| modcl D = 4 + — g-_ . (11 - 2x) . (114) 

So, the rule A + + A_ = 4 is broken in model D when k ~ 0(1). When k ^> 1, the equation 
can be written in the suggestive form 

A + +A -= 4+ (i^ (11 - 2 ^ (115) 
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where = 2M 2 is the number of scalars. We see that the violation of the rule A + + A_ = 4 
occurs when the number of scalars is of the same order as the number of color degrees of 
freedom, N 2 . 

Recall that the rule A + + A_ = 4 can be understood from AdS/CFT correspondence: A 
is related to the mass square m 2 of the bulk scalar by the equations A(A — d) = m 2 R 2 . How 
do we understand the fact that this rule is violated when there are 0(N 2 ) scalars? In fact, 
it is easy to come up with a mechanism leading to this effect within holography. Recall the 
AdS radius R is determined by the cosmo logical constant. Changing the boundary condition 
for the scalar field alters the vacuum energy (Casimir energy) associated with that scalar 
field (12] . The change in the vacuum energy is small for one scalar, but becomes of order 
one for 0(N 2 ) scalars. Thus we have 



d d 2 

A ± = -± ] /-+m 2 Rl (116) 

and, if R + /R_ = 1 + O^/N 2 ), then A + + A_ = 4 + O^/N 2 ), which is exactly what is 
found. We can construct an holographic model where the deviation of A + + A_ from 4 can 
be explicitly computed (see Appendix). Interestingly, in this simple model A + + A_ > 4, 
as in model D. 

The holographic model also lends support to the hypothesis that model C, with full 
chiral symmetry, exists in strong coupling. Indeed, the reason model C does not exist at 
weak coupling, from the holographic point of view, is that flipping the boundary condition 
for 2N 2 scalars from close to z 3 to close to z 1 is too much a disruption for the AdS geometry 
(for example, in terms of the change of the cosmological constant). However, when both A + 
and A_ are close to 2, the change of the vacuum energy is parametrically small in A + — A_ 



(see Eq. (All)), and flipping the boundary condition from z A+ to z A ~ is no longer a large 
disruption. Hence, the theory where all fermion bilinears have dimension A_, i.e., QCD*, 
should exist near the merger point. However, arguments based on holographic models can 
only be taken as suggestive at this moment. 



VII. CONCLUSION 

There have recently been several lattice studies seeking to find the boundaries of the 
conformal window in QCD |33 IH SSJ S§] and in other QCD-like theories 03 SSI US EDI ED 
E2J ESI El]- Interest in the phase transition between conformal and nonconformal theories 
is motivated in part by the invocation of approximate conformal symmetry in numerous 
theories for physics beyond the Standard Model. 

In this paper we have investigated the nature of such a phase transition, and we suggest 
that there is a wide class of theories where it is due to the merger and annihilation of fixed 
points. Several explicit examples of this phenomenon were given, and we speculate that this 
mechanism is also responsible for the chiral phase transition in QCD in the large N c , large 
Nf regime, at some critical value for Nf/N c . We show that this mechanism leads to the 
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BKT scaling behavior of the chiral condensate at the phase transition, and also implies the 
existence of the conformal theory QCD* which annihilates with QCD at the lower end of the 
conformal window. We tried, unsuccessfully, to construct QCD* in the perturbative regime, 
and argued that it should exist in the nonperturbative regime. It would be interesting to 
search for evidence for QCD* on the lattice. 

The models considered in the last section of our paper, in an attempt to find the UV 
fixed point of QCD, may be of interest in their own right. For example, these models may 
be used to explicitly realize the "unhiggs" [55], which behaves at high energies as a field 
with noninteger scaling dimension. 

The picture of the chiral phase transition realizes walking technicolor when Nj / N c is only 
slightly below the phase transition. In the holographic interpretation, conformality is lost 
when the mass squared of a bulk scalar drops below the BF bound. A naive application of 
AdS/CFT rules implies that the dimension of the operator ipif) is equal to 2 at the phase 
transition. This feature is explicit in the holographic model considered in the Appendix. This 
conclusion is, interestingly, in agreement with result from the Schwinger-Dyson approach, 
and also with Ref. [56J. The result illustrates that the dimension of the fermion bilinear on 
the IR stable branch cannot approach the unitarity bound used in Ref. [57] for estimating 
the maximal extension of the conformal window. 
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APPENDIX A: CASIMIR EFFECT IN ADS 5 

We consider a holographic model for a UV-IR pair of conformal theories. Each theory 
has a set of scalar operators that all have the scaling dimension 1 < A_ < 2 in one theory 
and 2 < A + < 3 in the other. Both theories are described by the same holographic dual 

n-V G -\^({dfaf + m 2 <%)) . (Al) 
i=i ' 

The two different theories correspond to the two choices of the boundary conditions on the 
scalar fields fa. The solution of the classical equations of motion is given by all fa = and 
AdS 5 metric: 

ds 2 = Rl z~ 2 (dz 2 + dx 2 - dt 2 ) , (A2) 

with 

Rl = -12 /V . (A3) 
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The loop expansion is controlled by dimensionless parameter k 2 Rq 3 ~ N 2 , where N c is the 
number of colors in the dual theory. One-loop contribution is not negligible in the N c — > oo 
limit if is also large, i.e., ~ N 2 . For simplicity we assume n^/N 2 <C 1 and compute 
A + + A_ — 4 to leading order in n^/N 2 . We use the technique of Ref. [42J. 7 The one-loop 
contribution of the scalar fields depends on the boundary condition. This contribution shifts 
the vacuum energy Vq to V±. Consequently, the AdSs curvature radius Rq becomes at one 
loop 

R 2 ± = -12/V±. (A4) 

It is convenient to measure lengths in units of R$, i.e., Rq — 1,Vq — —12, etc. 

The calculation of the vacuum energy is easier to perform after the Wick rotation t — > 
— 1x4. The correction to the vacuum energy is equal to 

J <fx E ^{V ± - V ) = 2k 2 V± logdet ± [(-V| + m 2 )^] , (A5) 

where index E denotes objects defined using the Wick rotated metric ds 2 E = (dz 2 +dx 2 +dx 2 ). 
The expression in the r.h.s. of Eq. (A6) is formal, since it is UV divergent. The derivative 
with respect to m 2 eliminates some but not all divergences: 

/d 
<f>x E ^—V ± = K 2 n 4> ti ± {[(-Vl + m 2 )^]" 1 ^}. (A6) 

The kernel of the operator [(— V|. + ^i 2 )^/(]e\ 1 is the Green's function defined by the 
following equation 

[-d z z- 3 d z + Q 2 z~ 3 + m 2 z~ 5 } G ±u {z, z'; Q) = 5{z - z') . (A7) 

where ±v refer to two different boundary conditions at z = 0: G± u ~ z 2±u , where v = 
\/4 + m 2 , and Q 2 = q 2 E . The solution to Eq. (A7) is 

G ±u (z, z'- Q) = z 2 z' 2 I ±v {Qz) K ±u (Qz') 9{z' - z) + (z <-> z') . (A8) 

Collecting results so far, we find 

^ v±=k2 ^ S^ g±v{z ' z]Q) - m 



The right hand side of Eq. (A9) is still UV divergent. However, the difference V + — V_ is 
finite: 



d 



K 2 n<p dm 2 



1 f°° 

(V+ - V.) = — / dQ Q 3 z 4 (I u (Qz) K u (Qz) - (u ~ -v)) 

1 2sinz/vr f°° 1 1 U) 

dxx 6 Kfix) = -——v(\ - ir) . 



8vr 2 vr L vy ' 12tt 



For earlier calculations of the Casimir energy in AdS space see Refs. [551 155] , 
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Using the fact that both boundary conditions are the same at v = 0, and thus V + 
v = 0, as well as dv 2 = dm 2 , we can write 



K 2 n^ 



(V + - VL) 



12tt 2 



dz> 2 v{\ - v 2 ) 



1 / V A V" 



6tt 2 V 3 



< 0. 



(Note that < v < 1.) 



Using Eqs. (A3), (A4) and in the regime \R± — R \ <C -Ro = b we can write 

1 



6 



(Vf-V-) 



The one-loop corrected scaling dimensions of the scalar operators become 



A± = 2± x U + m 2 R 2 ± . 



and thus 



U_ at 



(All) 



(A12) 



(A13) 



^.-4=^(^-^ = ^(4-^(1-^) >0. (A14 ) 

Since k 2 ~ A^ 2 the deviation of A + + A_ from 4 is O^^/N 2 ). At the point of merger 
v — > 0, A + + A_ = 4, hence A + = A_ = 2, with no correction of order 0(n<f)/N 2 ). 
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